clear, close, clc;
Nt = 50;
v = 1000;
h1 = 1000;
p1 = 50; p2 = 100;
u2 = p2 / p1;
i = sqrt(-1);


for j = 1:Nt
   W(j) = 10 ^ (j / 10);
   lanbta(j) =  2*pi*sqrt(2) / sqrt(W(j) /p1);
   A(j) = (tanh( sqrt(W(j) /p1)*h1 * exp(-i*pi/4) +atanh(sqrt(u2)))) ^ 2;
   
   B(j) = abs(A(j));
   C(j) = angle(A(j));
end
figure(1)
loglog(lanbta(:),B(:))
figure(2)
loglog(lanbta(:),C(:))